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0\ ; Abstract 

The Higgs boson may decay mainly to an invisible mode characterized by missing energy, instead 

(N 

J> ' of the Standard Model channels. This is a generic feature of many models where neutrino masses 

' 

■ arise from the spontaneous breaking of ungauged lepton number at relatively low scales, such 

I as spontaneously broken R-parity models. Taking these models as framework, we reanalyze this 

O ' 

' striking suggestion in view of the recent data on neutrino oscillations that indicate non-zero neutrino 



masses. We show that, despite the smallness of neutrino masses, the Higgs boson can decay mainly 

to the invisible Goldstone boson associated to the spontaneous breaking of lepton number. This 

D , requires a gauge singlet superfield coupling to the electroweak doublet Higgses, as in the Next 

^ ■ 

to Minimal Supersymmetric Standard Model (NMSSM) scenario for solving the /u-problem. The 
^ , search for invisibly decaying Higgs bosons should be taken into account in the planning of future 

' accelerators, such as the Large Hadron Collider and the Next Linear Collider. 
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I. INTRODUCTION 



Understanding the origin of mass is the main open puzzle in particle physics today. In 
the Standard Model all masses arise as a result of the spontaneous breaking of the SU(2) 
(E) U(l) gauge symmetry. This implies the existence of an elementary Higgs boson, not yet 
found. Stabilizing the mass of the Higgs most likely requires new physics and supersymmetry 
has so far been the leading contender. Another aspect of this problem is the smallness of 
neutrino masses. Despite the tremendous effort that has led to the discovery of neutrino 
mass 

HQS 

the mechanism of neutrino mass generation will remain open for years to come 
(a detailed analysis of the three-neutrino oscillation parameters can be found in m). The 
most popular mechanism to generate neutrino masses is the seesaw mechanism P,|6 J7,y,0|. 
Although the seesaw fits naturally in SO(IO) unification models, we currently have no clear 
hints that uniquely point towards any unification scheme. Therefore it may well be that 
neutrino masses arise from garden-variety physics having nothing to do with unification, 
such as certain seesaw vaHan. Q. a.d models w.h .adia^ve gene..atio„ QQ. In such 
models the physics of neutrino mass would then be characterized by much lower scales , 
potentially affecting the decay properties of the Higgs boson. This is especially so if neutrino 
masses arise due to the spontaneous violation of ungauged lepton number. In this broad class 
of models the Higgs boson will have an important decay channel into the singlet Goldstone 
boson (called majoron) associated to lepton number violation [l^ . 



h^JJ. 



(1) 



Here we focus on the specific case of low-energy supersymmetry with spontaneous viola- 
tion of R-parity, as the origin of neutrino mass. R-parity is defined as^j, = (_ i)3-B+l+25 
with 5*, B, L denoting spin, baryon and lepton numbers, respectively 15]- In this model 
R-parity violation takes place "a la Higgs" , i.e., spontaneously, due to non-zero sneutrino 
vacuum expectation values (vevs) jlfil. llTlllsl . In this case one of the neutral CP-odd scalars 
is identified with the majoron. In contrast with the seesaw majoron, ours is characterized 
by a small scale (TeV-like) and carries only one unit of lepton number. This scheme leads 
to the bilinear R-parity violation model, the simplest effective description of R-parity viola- 



tion [19] (for calculations including also trilinear terms see, for example 201 121 [ 



22 



23 



le model 



not only accounts for the observed pattern of neutrino masses and mixing 
but also makes predictions for the decay branching ratios of the lightest supersymmetric 



particle 
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from the current measurements of neutrino mixing angles 
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In previous studies j30| it was noted that the spontaneously broken R-parity (SBRP) 
model leads to the possibility of invisibly decaying Higgs bosons, provided there is an SU(2) 
® U(l) singlet superfield $ coupling to the electroweak doublet Higgses, the same that 
appears in the NMSSM. 

In this paper we reanalyse this issue taking into account the small masses indicated by 
current neutrino oscillation data Q] ^. We focus on the lowest-lying neutral CP-even scalar 
boson of the model. We show exphcitly that the presence of the SU(2) U(l) singlet 
superfield $ plays a triple role: (i) it gives a model where neutrino masses are obtained 
from first principles without any type of fine-tuning, even when radiative corrections are 
negligible, (ii) it solves the /i-problem "a la NMSSM" ^, and (iii) it makes the invisible Higgs 
boson decay in Eq. potentially the most important mode of Higgs boson decay. The 
latter is remarkable, given the smallness of neutrino masses required to fit current neutrino 
oscillation data. We also verify that the production of such Higgs boson in e~^e~ annihilation 
can be as large as that characterizing the standard case, and that therefore this situation 
should be taken as part of the agenda of future accelerators probing the mechanism of mass 
generation. 



II. MODEL WITH SPONTANEOUSLY BROKEN R PARITY 

The most general superpotential terms involving the Minimal Supersymmetric Standard 
Model (MSSM) superfields in the presence of the SU(2) ® U(l) singlet superfields (uf,Si, $) 
carrying a conserved lepton number assigned as (—1, 1,0), respectively is given as [3^ 

+ h'^^u^Sj + M'^u^S, + ^M^P + (2) 

The first three terms together with the fi term define the R-parity conserving MSSM, the 
terms in the last row only involve the SU(2) ® U(l) singlet superfields {Df,Si,^) ^, while 
the remaining terms couple the singlets to the MSSM fields. We stress the importance of 
the Dirac- Yukawa term which connects the right-handed neutrino superfields to the lepton 
doublet superfields, thus fixing lepton number. 

^ Ref. [sof assumed MeV-scale for the heaviest neutrino mass, inconsistent with the atmospheric data which 

points towards rrii, ~ 0.05 eV. 
^ Provided domain walls are either eliminated by imposing a Z2 -R-symmetry on the non-renormalizable 

operators , or that they are simply inflated away. 
^ The term linear in $ has been included in the first row as it is relevant in electroweak breaking. 
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A. Spontaneous Symmetry Breaking 



The presence of singlets in the model is essential in order to drive the spontaneous viola- 
tion of R parity and electroweak symmetries in a phenomenologically consistent way Like 
all other Yukawa couplings hu,hD,hE we assume that hi, is an arbitrary non-symmetric 
complex matrix in generation space. For technical simplicity we take the simplest case with 
just one pair of lepton-number-carrying SU(2) (8) U(l) singlet superfields, and S, in 
order to avoid inessential complication. This in turn implies, hij h and h]^ —>■ h^. 

The full scalar potential along neutral directions is given by 

Vtotai = \h^S + hiuiHu + MrS\^ + \ho^Hu + ilH^\^ + \h^v- + Mrv-\^ (3) 
+1 - ho^Ha - fiHa + Ku.uf + \ - hoH^Ha + hu^S -5'' + + ^^^\^ 

3 
i=l 

-Cs5'^ + Bm^Mru<^S + ^Sm,M$$2 ^ 1aaA$^ + /i.e.] 

3 

+ E + 1(9' + /) {\Hu\' - w - E l^^l')'' 

a 1=1 

where denotes any neutral scalar field in the theory. 

The pattern of spontaneous symmetry breaking of both electroweak and R parity sym- 
metries works in a very simple way. The spontaneous breaking of R parity is driven by 
nonzero vevs for the scalar neutrinos. The scale characterizing R parity breaking is set by 
the isosinglet vevs 

and 

m - ^. (5) 

We also have very small left-handed sneutrino vacuum expectation values 

fe>-^^ (6) 

The spontaneous breaking of R-parity also entails the spontaneous violation of total lepton 
number. This implies that one of the neutral CP-odd scalars, which we call majoron, and 
which is given by the imaginary part of 

f^Li / rr zj ^ .^^Li ~ .Vs Vr 

-^{vuH^ - vaHa) + 2^^^i+yS- y^^^ (7) 
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remains massless, as it is the Nambu-Goldstone boson associated to the breaking of lepton 
number. Note that this majoron is quite different from the one that emerges in the seesaw 
majoron model, as it is characterized by a different lepton number (one unit instead of two) 
and by a different scale, determined by the combination V = ^/f |. + f | ~ TeV . Note that 
Eq. is the origin of lepton number violation in this model and plays a crucial role in 
determining the neutrino masses. 

On the other hand, electroweak breaking is driven by the isodoublet vevs (Hu) = and 
(Hd) = with the combination = + v'^ + J^i'^li fixed by the W mass 



= (8) 



while the ratio of isodoublet vevs yields 

tan/3 = — . (9) 

This basically recovers the standard tree level spontaneous breaking of the electroweak 
symmetry in the MSSM js^ ^. 



B. Neutrino masses 



Since neutrino masses are so much smaller than all other fermion mass terms in the 
model, once can find the effective neutrino mass matrix in a seesaw-type approximation. 
From the full neutral fermion mass matrix, see Eq. ()A2jl . one calculates the effective 3x3 
neutrino mass matrix (m^^) as 

= -MSMh^Md, (10) 

where Mh is the 7x7 matrix of all other neutral fermion states, see Eq. ()A2j) . and the 3x7 
matrix m'^o is given as 

MS = ( mjo, mo ) , (11) 

where the matrices m^Oj, and mo are given in Eqs. ()A4|) and ()A7|) . The inverse of Mh is 
too long to be given explicitly here. 

After some algebraic manipulation, the effective neutrino mass matrix can be cast into a 
very simple form 

{v[vfj)ij = aAiAj + b{eiAj + ejAi) + ce^ey (12) 
^ We have verified explicitly, however, that radiative electroweak breaking may also occur. 
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where one can define the effective bilinear R-parity violating parameters and Aj as 
and 

K = GiVd + fJ^VLi (14) 

Here the parameter fi is 

^ = f^ + ho^, (15) 
while the coefficients appearing in Eq. ()12j) are given by 

(m^MRi-h^RVsfi + M^Mrii + hlMRVdVu)) (16) 



4/iDet(MH 

1 / 



8/iDet(MH) 



(hom^MR{hQMR + h^i)vu{vl - t;^)) (17) 



■ ((/loMij + /i/i)'i;^(2MiM2/i - m^VdVu)) (18) 



4/iDet(MH; 
and Det(MH) is given as 



Det(MH) = ^MR^miM2ii{M^MRii-h?iiVRVs + hlMRVdVu) 

- (^AfiVdimMR - h\RVs)vu + hlMRivj + vlf^ } (19) 

Note that Mr and M$ above are defined as 

Mr = Mr + h^, M$ = M$ + A^. (20) 
v 2 V 2 

The "photino" mass parameter is defined as = g^Mi + g''^M2. 

Eq. ()12p resembles very closely the corresponding expression for the explicit bilinear R- 



yery closely tne corre 



parity breaking model [ISl 120, |21|, |22, |2^ |2^ , once the dominant 1-loop corrections are taken 
into account. Note that the tree- level result of the explicit bilinear model can be recovered 
in the limit Mr, M$ oo. In this limit the coefficients b and c go to zero, while 

= 4Det(M^o) (^^^ 
In this limit only one non-zero neutrino mass remains. Whether the 1-loop corrections or 
the contribution from the singlet fields are more important in determining the neutrino 
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masses depends essentially on the relative size of the coefficient c in Eq. (fT^ compared to 
the corresponding 1-loop coefficient. Both extremes can be realized in our model. We note, 
however, that as discussed below large branching ratios of the Higgs into invisible final states 
require sizeable values of h and (as well as singlets not being too heavy). For such choices 
of parameters we have found that the "singlino" contribution to Eq. ()12|1 is usually much 
more important than the 1-loop corrections to the neutrino masses. 

Note also that the model does not predict whether the atmospheric (solar) mass scale is 
mainly due to the first (third) term in Eq. (jl2j) or vice versa. We have checked numerically 
that both possibilities can be realized and "good" points (in the sense of being appropriate 
for neutrino physics) can be found easily in either case. 



C. Scalar Mass Matrices 



With the above choices and definitions we can obtain the neutral scalar boson mass 
matrices as in Ref. by evaluating the second derivatives of the scalar potential in Eq. Q 
at the minimum. This results in 8 x 8 mass matrices for the real and imaginary parts of 
the neutral scalars ^. We have checked, in particular, that in the CP-odd sector we find 
both the Goldstone "eaten" by the Z° as well as the Goldstone boson corresponding to 
the spontaneous breaking of R-parity, namely the majoron, Eq. (0). In the basis A'q = 
{Hf, D^^, u'^^, z>^^, S^, D""^) these fields are given as. 

Go = {NoVd, -NoVu, NqVli, NoVl2, NoVl3,0,0,0) 

J = Ni{-N,Vd, NiVu, N2VL1, N2VL2, N2VL3, 0, Nsvs, -NsVr) (22) 

where the normalization constants Ni are given as 

1 



No 



Ni = vl^ + + vis 



N2 = vl + vl 
N3 = Ni + N2 



= (23) 



^ As already mentioned we assume, for technical simplicity that we have just one pair of lepton-number- 
carrying SU(2) U(l) singlet superfields. 
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and can easily be checked to be orthogonal, i. e. they satisfy Go ■ J = 0- 

In order to study the phenomenology of the scalar sector we need some information about 
the parameters of the SBRP model. Broadly speaking there are four types of parameters 
that are to a large extent undetermined. First there are Yukawa couplings, like h, ho and A. 
In contrast to hjj, hr, and He these are not fixed by fermion masses. Then there are MSSM 
parameters such as tan/5, the effective Higgsino mixing parameter /x, the supersymmetry 
breaking scalar mass parameters mo and Aq. These are partially restricted by negative 
collider searches for supersymmetric particles [34]. Then there are singlet sector mass pa- 
rameters, such as M/j, M$ and 5^. Finally there is the important Yukawa coupling hj^, which 
determines the strength of effective R-parity breaking parameters, through Eq. (|T3|l . This 
is constrained by neutrino oscillation data. In Section HYl we will discuss our strategy to 
choose the parameters in such a way that the results can be easily interpreted. We will also 
show there that a fully cubic superpotential, without any mass scale parameter such as the 
jlHuHd term, also leads to a realistic model consistent with neutrino oscillation data. 
Before that, however, we consider the corresponding Higgs boson phenomenology, focusing 
on Higgs boson production and decays, and stressing the potentially large invisible decay 
branching ratio. 



III. HIGGS BOSON PRODUCTION AND DECAYS 

Supersymmetric Higgs bosons can be produced at the e^e~ collider through their cou- 
plings to Z, via the so-called Bjorken process. In our SBRP model there are 8 neutral 
CP-even states Hi and 6 neutral CP-odd Higgs bosons Aj, in addition to the majoron J. 
One must diagonalize the scalar boson mass matrix in order to find the coupling of the 
massive scalars to the Z. The Lagrangean is 

8 

Chzz = Y^iV^Gpf'^MlZ^Z^^T^^H, (24) 

i=l 

with each rji given as a weighted combination of the five SU(2) (g) U(l) doublet scalars, 

3=1 

where is the 8x8 rotation matrix for the CP-even scalars. Note that we leave the discus- 
sion of the CP-odd scalars for elsewhere. Moreover, here we focus mainly on the production 
of the lightest CP-even supersymmetric Higgs boson h = Hi. The main difference between 
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the production of this state and the hghtest CP-even Higgs boson of the MSSM is the fact 
that ours contains an admixture of the SU(2) (8) U(l) singlet scalar fields z/"^ and S, and its 
coupling to the Z is correspondingly reduced by a factor 



?7 = ?7i < 1 



(26) 



in comparison with the Standard Model case ^. When the lightest CP-even Higgs boson is 
mainly singlet its production cross section in e'^e" annihilation will be suppressed. 

We now turn to the lightest Higgs boson decays. Given that other MSSM decay modes 
are less important, we are particularly interested here in the ratio 

T{h J J) 



R. 



Jb 



(27) 



r(h hh) 

of the invisible decay to the Standard Model decay into b-jets. For this we have to look 
separately at the decay widths, 



T{h J J) 



9hjj 
327rmh 



and 



r{h hh) 



1-4 



rrih 



3/2 



(28) 



(29) 



Stt cos^ (5 

Prom these expressions we see that r(/i — > hh) will be small if the component of the lightest 
Higgs boson along is small. On the other hand the magnitude olV{h ^ J J) will depend 
on the Qhjj coupling. This is in general given by a complicated expression, but for the 
situation that we are considering here with 



VLi < Vd, Vu < Vr, Vs 



we have to a very good approximation 



J -(0,0,0,0,0,0,^,-^) 



(30) 



(31) 



where = f | + vj^. Under this approximation we can write the coupling g'- for the vertex 
h[JJ of the Majoron with the unrotated Higgs boson h'^, in the following form 



vsVr 



For the MSSM we have a reduction given by 77 = ^R'li + ^Ri2 = sin(/? — a) in the usual notation. 
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' 2e,_2^^^^^^ (z = 3,4,5) 



y2 . 



96 



9r - 

9's - -h^P (32) 

where Mr and M$ have been defined in Eq. ()2()j) . 

From these expressions we conclude that Qhjj can be large in two situations. The first 
is, of course, if the lightest Higgs boson is mainly a combination of the and 5* fields. In 
this case not only g^jj will be large, but also r(/i — > bb) will be small suppressing h — > bb. 
Unfortunately the production would be suppressed, as singlets do not couple to the Z. 
The phenomenologically novel and interesting situation is when h and ho are large. In this 
case the Higgs boson behaves as the lightest MSSM Higgs boson (with moderately reduced 
production cross section) but with a large branching to the invisible channel h —>■ J J. 

The sensitivities of LEP experiments to the invisible channel h —>■ J J have been discussed 



since long ago [36|, |37| and the current status has been presented in Ref. [38|. In order to 



evaluate the experimental sensitivities to the parameters of the model we must take into 
account both the production as well as Higgs decays. 



IV. NUMERICAL RESULTS 



In this section we discuss the numerical results on the invisible decay of the Higgs boson 
in our model. We start with a brief discussion of the SBRP parameters. 

Unknown parameters of the spontaneous R-parity breaking model fall into three different 
groups. First, there are the MSSM parameters, mainly the unknown soft SUSY breaking 
terms. The second group of parameters are the and left-handed sneutrino vevs vl,- We 
trade the latter for the parameters Aj using Eq. These six parameters occur also in the 
explicit bilinear model. And, finally, there are the parameters of the singlet sector, namely 
singlet vevs vr, vs and Yukawa couplings h, ho and A and the singlet mass terms Mr, 
M$, 5^, as well as the corresponding soft terms. 
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We have checked by a rather generous scan that the results presented below qualitatively 
do not depend on the choice of MSSM parameters, as expected. Thus, for definiteness we 
will fix the MSSM parameters in the following to the SPSla benchmark point j^^, defined 

by 

mo = 100GeV mi/2 = 250GeV tan/5=10 Aq = -lOOGeV fi< (33) 

We have run down this set of parameters to the electro-weak scale using the program package 
SPheno We stress again that different choices of MSSM parameters will not lead to 



qualitatively different results. 



A. General case 



We first consider the general model defined by the superpotential in Eq. Q reduced to 
one generation of z/'^ and S fields. For the singlet parameters we choose as a starting point 
vj^ = = = -150 GeV and Mr = -M$ = 5 = 10^ GeV, as well as h = 0.8, ho = -0.15 
and A = 0.1. We have tried other values of parameters and obtained qualitatively similar 
results to the ones discussed below. 

The explicit bilinear parameters are then fixed approximately such that neutrino masses 

f] nnn 

and mixing angles |^ are in agreement with experimental data [H, 12, |3||. Slightly different 
values of parameters are found, depending on whether the first or the third term in Eq. (|12p 
is responsible for the atmospheric neutrino mass scale. Both possibilities lead to very similar 
results for the invisible decay of the Higgs. This can be understood quite easily. The ratio 
of the atmospheric and solar neutrino mass scale is only of the order of (4 — 7) ^ and the 
changes in parameters A and e are only of the order of the square root of this number. 
Such a small change can always be compensated by a slight adjustment of other parameters, 
leading to the same (or very similar) final result. 

After having defined our "preferred" choice of parameters in the following we will vary 
one unknown parameter at a time. We now turn to a discussion of the results. In Fig. (Q) we 
show the ratio Rjh as a function of ■i]'^ for different choices of h (left) and for different choices 
of Vr (right) and all other parameters fixed. Larger values of Rjb are found for smaller 
values of t], as expected. However, one sees explicitly that even for values of ?7 ~ 1, Rjb can 

^ In a hierarchical model, such as the one discussed here, the square roots of the Amfj are approximately 
equal to the larger mass. 
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be larger than 1. This means that the hghtest Higgs can decay mainly invisibly, even when 
the cross section for its production is essentially equal to the usual (MSSM) doublet Higgs 
boson cross section. This is the main result of this work. 
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Figure 1: Ratio Rjh-, defined in Eq. as function of rf. a) to the left, for different values of 

the parameter /i, from top to bottom: h = 1, 0.9, 0.7, 0.5, 0.3, 0.1. b) to the right, for different 
values of the parameter vr = vs: -vr = 150, 200, 300, 400, 600, 800, 1000 GeV. The plots show 
explicitly that Rjb > 1 is possible even for 77 ~ 1. This is the main result of the current paper. 

In Fig. (j21) we show Rjb as function of = + Vg (left) and as function of h (to 
the right). The figure shows that large values of Rjb are obtained for small values of V 
and for large values of h. The decreasing of Rjb with increasing values of V can be easily 
understood, since in the limit V —>■ 00 the majoron should obviously decouple. 
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Figure 2: Ratio Rjh, defined in Eq. (j2Z|i as function of V (left) and as function of h (to the right). 
Small (large) values of V {h) lead to large values of Rjb. 



Other singlet-sector parameters also can have an important impact on Rji,, as demon- 
strated in Fig. 0. As shown in the left panel of this figure, larger values of ho lead to larger 
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values of Rjb- For values of h smaller than about h ^ 0.75 (for our specific choice of the 
other parameters) the order of the lines is exchanged. This is due to a level-crossing in the 
eigenvalues. Below this value, the lightest Higgs is mainly a singlet and thus even though it 
decays dominantly invisibly its production cross section is very much reduced. 

On the other hand, the right panel of Figure (jH]) shows that the value of f$ is normally 
somewhat less important than the value of V in determining Rjb. Again this can be quali- 
tatively understood since V is the parameter whose magnitude determines the breaking of 
lepton number (indeed, with the help of the approximate couplings g'^ in Eq. (j32|) one can 
see that the parameters h, h^, vr and vs should be the most important ones). 




Figure 3: Ratio Rjb-, defined in Eq. as a) left figure: function of \h\ for —h^ = 

0.3, 0.1, 0.03, 0.01, 0.001 (on the right part of the plot from top to bottom). The right panel 
b) gives Rjb as function of for different values of the parameter vr = vs for —vr = 
150, 175, 200, 300, 400, 600, 800, 1000 GeV. 

As a summary of this section we conclude that large branching ratios of the doublet-like 
Higgs boson into invisible final states are possible in the SBRP model, despite the smallness 
of the neutrino masses indicated by oscillation data. Large values of Rjb occur for large 
values of the Yukawa couplings and for small values of vr. The presence of the field $ plays 
a crucial role in getting the invisible Higgs boson decays that are not suppressed by the 
small neutrino masses. 



B. Cubic only superpotential 

Before concluding we illustrate the results we have obtained for the case of a restricted 
SBRP .odel de^cHbed by tbe =upe.poten«aI .n E,. O coMa,n.g only cuMc te... Q. 
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The restricted model provides a potential "solution" to the /i problem in the context of 
spontaneous R-parity violation. We give results for the same parameter choices as above, 
except that no mass parameters are now present in the basic superpotential. 
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Figure 4: Ratio Rjb, defined in Eq. ()27() . as function of r/^, a) to the left, for different values of the 
parameter h and b) to the right, for different values of the parameter vr = vs- As in the general 
case (Fig. large values of Rjb can be found even for r/ ~ 1 also in the cubic-only case. 
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Figure 5: Ratio Rjb, defined in Eq. as function of the parameter V (left) and as function of 
h (right). The qualitative behaviour is similar to the general case, compare to figure |2j 

Even though acceptable physical solutions consistent with experiment (supersymmetric 
particle searches as well as neutrino oscillation data) are somewhat harder to find, they 
exist. Figs. (jH) and © show Rjb as function of r/^ and as function of h and V for the cubic- 
only case, compare to Figs. ((H) and Q for the general case. As can be seen, the qualitative 
behaviour is very similar in all cases, although the parameters for which acceptable solutions 
are found are usually restricted to narrower ranges in the cubic-only case. These figures 
demonstrate that also in the cubic-only case large production cross section and large invisible 
branching ratios for the lightest Higgs decay can occur at the same time. 
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V. DISCUSSION 



We have discussed the possibihty of an invisibly decaying Higgs boson in the context of 
the spontaneously broken R-parity model. One of the neutral CP-odd scalars in this model 
corresponds to the SU(2) ® U(l) singlet Nambu-Goldstone boson associated to the breaking 
of lepton number. In contrast to the MSSM, where the Higgs boson can decay invisibly only 
to supersymmetric states (in regions of parameters where the Higgs is heavier than twice 
the lightest neutralino mass) in our case the Higgs can decay mainly due to Eq. instead 
of the Standard Model channels, over large regions of parameters, given that there is no 
kinematical barrier for this decay. We have reanalysed this striking suggestion in view of 
the recent data on neutrino oscillations that indicate non-zero but small neutrino masses. 
We have explicitly shown that (i) despite the smallness of neutrino masses, invisible Higgs 
boson decay may indeed provide the most important mode of Higgs boson decays and (ii) 
its production cross section need not be suppressed with respect to that characterizing the 
standard MSSM case. As a result, our analysis indicates that invisibly decaying Higgs 
bosons should be an important topic in the agenda of future accelerators, such as the Large 
Hadron Collider and the Next Linear Collider. In fact the interest on this possibility goes 
beyond the model we have taken as framework, it is much more general. However the SBRP 
model provides an attractive explanation for the origin of the neutrino masses that can 
also be probed at future collider experiments through the predicted pattern of the Lightest 
Supersymmetric Particle (LSP) decays which directly traces the experimentally observed 
neutrino mixing angles. 
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Appendix A: NEUTRINO-NEUTRALINO-SINGLINO MASS MATRIX 



In the basis 

{-iX', -iX^, Hd, Hu, z/e, i^f,, i^T, S, $) 
the mass matrix of the neutral fermions following from Eq. Q can be written as 



(Al) 



/ M^o m^o^ m^Oj,c m^o^ \ 



M 



N 



mjo^ mo 







M^g Ms^ 



(A2) 



where the matrix M^o is the MSSM neutralino mass matrix: 







(A3) 



Ml 

M2 +lgvd -\gvu 

- la'vd +lgvd -fx 

\ + \g'vu -\gvu -fi 

Here, jj, = fi + hQV^/\/2. m^oj, is the R-parity violating neutrino-neutralino mixing part, 
which also appears in explicit bilinear R-parity breaking models: 



( - \g'vLe \gVLe ^ 



T 



^g'vLpi \gvLi, 



(A4) 



y - \g'vLT \gvLT o y 

where vii are the vevs of the left-sneutrinos, are defined by = -^hlvR, and vr is the 



vev of the right-sneutrino. 
Here m^^oj^c is given as 



0, 0, 0, 



V2 



VLi 



(A5) 
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and m^o^ is 



m!^o^ = (0, 0, --^hoVu, --^hoVd) (A6) 



The "Dirac" mass matrix is defined in the usual way 

1 

7! 



mD)i = -^Kvu (A7) 



The and S states are coupled by 
Mjc^ and Mj^ are 



(M,cs) = Mji + h^ (A8) 



Mjc^ = {{vs)) (A9) 

Ml* = (M) (AlO) 
Here, (f^j) = hvR and (^5) = hvs- Finally M* is 

M* = A'U + (All) 
v2 

We briefly comment on the case of three generations of neutral fermions in the singlet 
sector. For three copies of u'^ and S fields the mass matrix of the neutral fermions can be 
written in exactly the same form as given in Eq. ()A2|) with some rather straight-forward 
generalizations of the above definitions. These changes are: h and h]^ become 3x3 matrices 
h^^ and h^J . In Eq. f)A5|) the matrix becomes a 3 x 4 matrix, M/? is a symmetric 3x3 matrix 
and Eqs. ()A9|) and ()A10|) have to be replaced by 

M^.^ = i{vs,),{vs,),{vs,)) (A12) 

Mj^ = {{vm),{vR,)AvR,)) (A13) 

where (vm) = Y.j ^^'^f and {vs,) = Y.j h'^v^- 

Notice that even with three generations of z/^ and S fields one neutrino mass is zero at 
the tree-level. 
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Appendix B: THE NEUTRAL SCALAR MASS MATRIX 



The 8x8 scalar mass matrix is a symmetric matrix that in the basis of the real part of 
H^, z/j, $, S, u'^) can be written in the form, 



Mf-n Mf, Mf, 



(Bl) 



where M^j^ is a symmetric 2x2 matrix, M~~ and Mgg are symmetric 3x3 matrices, while 



and are 2x3 matrices and finally M~^^ is (a non-symmetric) 3x3 matrix. In 
this notation L denotes the sneutrinos and 5* the singlet fields. 

We can write the mass matrix by giving the components of the various blocks. We get. 



HH 



3 



i=l 

1 

Mi'^,, = ^(/ + ^>^ + f^cot/?-^^X^A.»,.-^M^^X^/.>,. (B4) 



^HH,, = --{g^ + g'^)vdVu-n + hlvuVd (B3) 



u 



i=l 

where, 

A 1 ^/2, 

Q = Bfi- 6^ho + -hovl + -hh^VRVs + —AhahoVi^, + —hoM^v.^, (B5) 



and /i, Mji and M$ are defined in Eqs. (fl3|) and ((201) ■ 



LL 



^> - ^ — K 2^ K^Lk - -i^Mr hi (B6) 
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1 1 



= ^MnVuK (B8) 



k=l 



Kh. = lla' + a^soLi-^fiVEhi (BIO) 

1 /2 /2'--- ^ 

^HL,i = -^l/ + + — + —MRVsK + VuhiJ2htvLk (Bll) 

fc=i 



<5n 



V2hoiJ,Vd - (Ah, + Vu - ^hoVR ^ h'^.VLk (B12) 

ifc=i 



= —hhoVRVu (B13) 
--hho vsVu - —H ^ h^VLk (B14) 



k=l 



-j^ 3 

= V2^oA*^^« - ^/io [Ah, + ^;<i + -/i ^;5 XI /i^^^Lfe (B15) 

k=\ 

mSs,, = -^/^/^o VRVa + 5^ /i^fc (B16) 



2 

fe=i 



1 ^ /2 ^ 

= - 9^^0 ^^5^^<i + v^vr X /l^/l^ + ^ ^/^>^^Lfc (B17) 



2 

fe=l A;=l 
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»Mfs 



Mis,, = 



Mi 



SS12 ' 



M: 



SS22 



Mil, 



where 



f$ 4 



— 7rh [Ah + M$) + —ho [Ah^ + M$) + -/iq 

2 2 1)$ 2 



ik=l 



fe=l 



V2 
2 

72 
2 



^ 1 ^ 

h (Ah + vr + \/2/i Mr^;5 +^hVuJ2 ^t^Lk 

k=l 

h (Ah + V5 - -hoVd hlvLk + V2h MrVr 



-r Mr > h^VLk 



Mls,, = T + h^VRVs 



Vr 2 Vr 



k=l 



VLk 



k=l 



1 1 a/2 

r = BmuMr - S'^h + -hXvl - -hhoVuVd + —h {Ah + M$) v^ 



(B18) 

(B19) 

(B20) 

(B21) 
(B22) 
(B23) 

(B24) 



Appendix C: THE NEUTRAL PSEUDO-SCALAR MASS MATRIX 

The 8x8 pseudoscalar mass matrix is a symmetric matrix that can be written in the 
form, 













Mi' 






LS 


M§7 




Mis 



(CI) 



where the blocks have the same structure as before. We can write the mass matrix by giving 
the components of the various blocks. We get, 
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M^H,^ = tan /? + ^/i^ 5^ /^t 

2 Vd 



« = 1 fc = l 

where Vt and /i are given in Eqs. ()B5|) and psp. 



2 t^Li " 2 VLi 



2 ^Li " " 2 VLi " 



k=l 



2 
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Mhs,, = -Y^o [^ho - M^) Vu + ^hoVR J2 hlvLk (CIO) 

k=l 

Mhs,2 = -\hhoVRVu (Cll) 

fc=i 

^£.1 = (^'^0 -M^)vd + -hvsY, ^tvLk (C13) 

fe=i 

^IL = - 2^^o ^^^-^ + 5Z ^'""Lk (C14) 

fe=i 

^£23 = -l^ho VsVd -^Yl "^hM^Lk (C15) 



2 4 

— Trh [Ah + M$) h — /lo (^/lo + ^ 0^0 > K^Lk 

2 2 2 ^ — ' 

k=l 

+26' X + Xho vuv, - Xh VRVs -Ih^f] htvLk - ^hMn^±^ (C16) 

2 ^-^ 2 

k=l 

3 



-^h (a, - VR - ^hv^Y.htvLk (C17) 

k=l 

-^h [a^ - vs - -hov, J2 htvLk (C18) 

k=i 

-r'^-^MZ-thtv^k (C19) 



k=l 

-r (C20) 
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k=l k=l 

where F is given in Eq. ()B24|) . 



[1] Super-Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 (1998), 
|liep-ex/ 9807003l. 

[2] SNO Collaboration, Q. R. Ahmad et al, Phys. Rev. Lett. 89, 011301 (2002), 
|mlcl-ex/020 4008l. 

[3] T. Araki et al., [KamLAND Collaboration], arXiv:hep-ex/0406035 Talk by G. Gratta at 
Neutrino 2004, 14 19 June 2004, Paris, |http://neutrino2004.in2p3.fr7 K. Eguchi et al, Phys. 
Rev. Lett. 90, 021802 (2003), |hep-ex/0212"02T| . 



[4] M. Mahoni, T. Schwetz, M. A. Tortola and J. W. F. Vahe, New J. Phys. 6, 122 



(2004), hep-ph/0405172 . See also M. Maltoni et al, Phys. Rev. D68, 113010 (2003), 
|hep-ph /0309130 , and references therein. 
[5] M. Gell-Mann, P. Ramond and R. Slansky, (1979), Print-80-0576 (CERN). 
[6] T. Yanagida, In Proceedings of the Workshop on the Baryon Number of the Universe and 

Unified Theories, Tsukuha, Japan, 13-14 Feb 1979. Edited by O. Sawada and A. Sugamoto. 
[7] R. N. Mohapatra and G. Senjanovic, Phys. Rev. D23, 165 (1981). 
[8] Y. Chikashige, R. N. Mohapatra and R. D. Peccei, Phys. Lett. B98, 265 (1981). 
[9] J. Schechter and J. W. F. Valle, Phys. Rev. D22, 2227 (1980). Phys. Rev. D25, 774 (1982). 
[10] R. N. Mohapatra and J. W. F. Valle, Phys. Rev. D34, 1642 (1986). 
[11] A. Zee, Phys. Lett. B93, 389 (1980). 
[12] K. S. Babu, Phys. Lett. B203, 132 (1988). 
[13] J. W. F. Valle, Prog. Part. Nucl. Phys. 26, 91 (1991). 
[14] A. S. Joshipura and J. W. F. Valle, Nucl. Phys. B397, 105 (1993). 

[15] C. S. Aulakh and R. N. Mohapatra, Phys. Lett. B 119, 136 (1982). G. G. Ross and 
J. W. F. Valle, Phys. Lett. B 151, 375 (1985). J. R. Ellis et al, Phys. Lett. B 150, 142 
(1985). L. J. Hah and M. Suzuki, Nucl. Phys. B 231, 419 (1984). 

[16] A. Masiero and J. W. F. Valle, Phys. Lett. B251, 273 (1990). 

[17] J. C. Romao, C. A. Santos and J. W. F. Valle, Phys. Lett. B288, 311 (1992). 

[18] M. Shiraishi, I. Umemura and K. Yamamoto, Phys. Lett. B313, 89 (1993). 

[19] M. A. Diaz, J. C. Romao and J. W. F. Valle, Nucl. Phys. B524, 23 (1998), hep-ph/9706315. . 

23 



E. J. Chun and S. K. Kang, Phys. Rev. D61, 075012 (2000), |hep-p h/9909'429l. 

A. Abada, S. Davidson and M. Losada, Phys. Rev. D65, 075010 (2002), |hep-ph/0111332] . 

M. A. Diaz et al, Phys. Rev. D68, 013009 (2003), hep-ph/030 2021. . 

M. Hirsch et al, Phys. Rev. D62, 113008 (2000), |hep-ph/0004Tl5l , Err-ibid. 
D65:119901,2002. 

J. C. Romao et al, Phys. Rev. D61, 071703 (2000), |hep-ph/9907"499] . 
M. Hirsch and J. W. F. Vahe, New J. Phys. 6 76 (2004), |hep-ph/0405015| 
M. Hirsch et al, Phys. Rev. D66, 095006 (2002), |hep-p h/0207334l. 



W. Porod et al, Phys. Rev. D63, 115004 (2001), |hep-ph/0011248| . 

D. Restrepo, W. Porod and J. W. F. Valle, Phys. Rev. D64, 055011 (2001), |hep-ph/0104"040] . 
M. Hirsch and W. Porod, Phys. Rev. D68, 115007 (2003), |hep^ph/0 307364 . 
F. de Campos et al, Nucl. Phys. B451, 3 (1995), lhep^ph/9502237 . 

C. Panagiotakopoulos and K. Tamvakis, Phys. Lett. B446, 224 (1999), hep-ph/ 9809475] . 
J. C. Romao, F. de Campos and J. W. F. Valle, Phys. Lett. B292, 329 (1992), 
|hep-ph/9207269| . 

R. Barbieri, S. Ferrara and C. A. Savoy, Phys. Lett. B119, 343 (1982). 

DELPHI Collaboration, P. Abreu et al, Phys. Lett. B502, 24 (2001), |hep-e x/0102045|. 

J. C. Romao, A. loannisian and J. W. F. Valle, Phys. Rev. D 55 (1997) 427 

|arXiv:hep-ph /960740T1 . 



A. Lopez-Fernandez et al, Phys. Lett. B312, 240 (1993), |hei>ph/9304255, . 
F. de Campos et al, Phys. Rev. D55, 1316 (1997), hep-ph/9601269 |. 

J. Abdallah et al [DELPHI Collaboration], Eur. Phys. J. C 32 (2004) 475 
|arXiv:hep-ex/0401022| . 

B. C. Allanach et al, Eur. Phys. J. C25, 113 (2002), |hep-ph/0202233l . 

W. Porod, Comput. Phys. Commun. 153 (2003) 275 |arXiv:hep-ph/0301101| 



24 



